Spontaneous formation of complex micelles from homogeneous solution 
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We present an extensive computer simulation study of structure formation in amphiphilic block 
copolymer solutions after a quench from a homogeneous state. By using a mesoscopic field-based 
simulation method, we are able to access time scales in the range of a second. A "phase diagram" 
of final structures is mapped out as a function of the concentration and solvent-philicity of the 
copolymers. A rich spectrum of structures is observed, ranging from spherical and rodlike micelles 
and vesicles to toroidal and net-cage micelles. The dynamical pathways leading to these structures 
are analyzed in detail, and possible ways to control the structures are discussed briefly. 
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Amphiphilic molecules in solution such as lipids or am- 
phiphilic block copolymers self-assemble into a variety of 
structures, e. g., spherical or cylindrical micelles, lamel- 
lae, and vesicles[l|. In the case of single-component am- 
phiphiles, vesicular and toroidal structures are energet- 
ically less favorable than lamellar and cylindrical struc- 
tures due to the energy penalty for bending. Neverthe- 
less, they may be stabilized by entropic or kinetic factors 
and still form spontaneously. These structures provide 
new opportunities for designing soft materials with en- 
hanced functionalities for various applications, such as 
complex micro release systems or templates for nanode- 
vice fabrication@,S,|i,[3,S,0|. A detailed understanding 
of the aggregation process is crucial to understand and 
eventually control their formation. 

In the past decades, a number of experimental stud- 
ies have revealed the rich diversity of micellar morpholo- 
gies displayed by amphiphilic systems. Besides spherical, 
rodlike or wormlike micelles, micelles with various special 
topologies have been observed, such as unilamellar and 
multilamellar vesicles, onion vesicles, genus vesicles 0,0, 
0, HI, 0, 0, S, B [131; toroidal micelles with one or sev- 
eral rings, and net or cage micelles 11 , 13, HI, 0]. The 
experiments indicated that the micellar structures and 
size distributions not only depend on molecular parame- 
ters, i.e., the chain length of the amphiphilic molecules, 
the hydrophilic-to-hydrophobic ratio, the molecular stiff- 
ness and the intermolecular interactions, but also on 
system parameters such as the concentration, and on 
kinetic factors such as the diffusion ability of the am- 
phiphilic molecules and the details of the manufacturing 
process. Unfortunately, detailed dynamical information 
on the process of spontaneous micelle formation is scarce. 
Only a few groups have captured the process of sponta- 
neous vesicle formation in solutions of amphiphile mix- 
tures and proposed a possible pathway of vesicle forma- 
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tion 

According to this "standard" pathway (PC), the am- 
phiphilic molecules first self-assemble into small spheri- 
cal micelles, these then coalesce to rods, the rods trans- 
form themselves to bilayers, and finally, the bilayers bend 



around and close up to vesicles. The last two steps are 
driven by the rim energy of the bilayers. The mechanism 
has been confirmed by computer simulations of different 
coarse-grained models [13, ED, H2, HI, 0, HI- It clearly 
contributes to the formation of vesicles in amphiphilic 
systems. However, it cannot explain the existence of 
complex toroidal structures, since no force pushes rod- 
like micelles with two detached end caps to form rings, 
let alone cage structures. In a recent paper, we have re- 
ported the existence of an alternative pathway of vesicle 
formation in copolymer solutions[26] (PG). In this path- 
way, the micelles do not coalesce, but simply grow by 
attracting copolymers from the solution. Once a critical 
micelle size is exceeded, copolymers start to flipflop such 
that the micelle core becomes solvent-philic ("semivesi- 
cle" state). Finally, solvent diffuses inside the core, and 
the semivesicle swells into a vesicle. The two pathways 
are illustrated in Fig. Q] (see below for simulation details). 
The pathway PG may provide a possible route to toroidal 
structures. 

In this letter, we report on an extensive systematic 
study of structure formation in a single-component am- 
phiphilic diblock copolymer system. The final self- 
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FIG. 1: Pathways of spontaneous vesicle formation in copoly- 
mer solutions. (a) PC: Micelle coalescence (al,a2), bi- 
layer formation (a3), and bending (a4); (b) PG: Micelle 
growth (bl,b2), internal reorganization to semivesicle (b3), 
and swelling into vesicle (b4). The parameters of the simu- 
lations are \bs = 0.128 and $ p = 0.2, (a) $p = 0.15 (b). 
Times t are given in units of T = 10 4 ro. Here and through- 
out the paper, structure snapshots show isodensity surfaces 
of ^-blocks at $a = 0.625. 



assembled structures depend strongly on the copolymer 
concentration and the interaction parameters. By vary- 
ing the latter over a wide range, we can map out the 
final topologies in a unifying "phase diagram". The sim- 
ulations allow to investigate the formation process in de- 
tail. Both pathways PC and PG described above can 
be observed, depending on the copolymer concentration. 
We find that vesicles and rodlike micelles may form ir- 
respective of the pathway, but toroidal structures only 
form via the pathway PG. Our results thus demonstrate 
that complex structure formation is not only controlled 
by the molecular packing parameters, but also, crucially, 
by the details of the segregation kinetics. 

Complex vesicle formation has also been studied re- 
cently by Sevink and ZvelindovskydlEzI]- They consid- 
ered copolymers made of two incompatible blocks that 
were both basically solvent-phobic. As a result, the 
copolymers aggregated to compact, internally structured 
droplets (onion vesicles) with a relatively low solvent con- 
tent (a few percent|27]). In contrast, in this work, we 
focus on copolymers with strongly solvent-philic compo- 
nents, and on open and hollow structures. 

We consider a system of amphiphilic diblock copoly- 
mers P (copolymer volume fraction <&p) with solvent- 
phobic blocks A (chain fraction ca) and solvent-philic 
blocks B (chain fraction cp = 1 — ca), immersed in a sol- 
vent sjHI]. The monomer interactions are characterized 
in terms of Flory-Huggins parameters xab, Xas, and 
Xbs- A compressibility modulus kh ensures that the lo- 
cal density (polymer plus solvent) is roughly constant. 
The time evolution of the system is modeled with Exter- 
nal Potential Dynamics a dynamic density functional 
theory which locally conserves densities and is approxi- 
mately valid for Rouse-type chain dynamics, but neglects 
hydrodynamics and reptation (see Refliil for a compila- 
tion of the dynamical equations). The relevant dynam- 
ical model parameters are the mobility coefficients D$ 
and Dp of the solvent and the copolymer. 

The parameter xbs and $p were variable. The other 
model parameters were set to xab = 0.896, xas = 1-024, 
K H = 1.176, cb/ca = 0.133, and D S /D P = 17. The high 
values of xas and xab ensured that the ^4-blocks segre- 
gate well from the S-blocks and the solvent. Copoly- 
mers with short -B-blocks were used to stabilize bilayer 
structures in the regime where B is swollen with solvent. 
The remaining parameters, Dp and R g (the unperturbed 
radius of gyration of the chains), set the length scale 
(r := Rg/3) and the time scale (r := r^/Dp) of the 
simulation. Mapping these to real copolymer solutions 
such as, e.g., those studied in Ref.l3(il (Dp ks 10~ 6 cm 2 s -1 
and R g w 30nm) , we find that our time unit corresponds 
to roughly to ~ 10~ 6 s. 

On the technical side, the parameters of the simula- 
tion were as follows: The time step for integration of the 
dynamical equation was chosen At = 0.02to, on a spa- 
tial grid with grid size rg = R g /3. The contours of the 
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FIG. 2: "Phase diagram" of final structures after a sudden 
quench from an initially homogeneous copolymer solution, for 
a range of solvent-philicities xbs and copolymer volume frac- 
tions $p. The final structures in the regions V1/V2 corre- 
spond to vesicles (f), RS1/RS2 to mixtures of rod and sphere 
micelles (d,e), Tl to ring micelles(c), and T2 to toroidal mi- 
celles (a,b). In the region R0, the solution stayed homoge- 
neous. The dotted lines are guides for the eye. The dashed 
line separates two dynamical regimes where the structure for- 
mation proceeds along different pathways: Micelle coalescence 
(pathway PC) in the regions RSI and VI, micelle growth 
(pathway PG) in the regions RS2, V2, Tl, and T2. The solid 
line shows the function x* = 1-3 — 11. $p (see text and Fig. [4] 
for explanation). 



copolymer chains were discretized with A = 17 steps. 
A small Gaussian noise was added to mimick the effect 
of thermal fluctuations [ifj. The longest total simulation 
time was 7.2 x 10 5 To, corresponding to 0.72 second. 

Fig. [2] shows a diagram of final structures, obtained 
after quenching the system suddenly from an initially 
perfectly homogeneous state. The solvent-philicity xbs 
ranges from positive to negative, in order to represent a 
wide class of amphiphilic block copolymers from nonionic 
to ionic. In addition, the copolymer volume fraction $p 
was varied in the region where interesting structures were 
observed. 

In order for structure formation to take place, the 
copolymer concentration $p must exceed a certain criti- 
cal value, which can be identified with the CMC (crit- 
ical micelle concentration). At lower copolymer con- 
centrations, the translational entropy of the copolymers 
prevents them from aggregating. The shape of the fi- 
nal structures depends on the solvent-philicity of the B- 
block, (—Xbs)- At moderate xbs (xbs > 0), the system 
favors bilayer structures and the copolymers aggregate to 
vesicles. As the solvent quality for the -B-block increases 
(xbs < 0), the B-block swells. Its radius Rb roughly 
scales with[ai| R B ~ (cpA0 3/5 (V 2 - xbs) 1/5 (N is the 
chain length) . This in turn decreases the critical packing 
parameter[l|, the bilayers become unstable and give way 
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FIG. 3: Formation of toroidal structures at \bs = —0.512: 
(a) Rings ($p = 0.2), (b) Toroidal micelles ($p = 0.18), (c) 
Cage micelles (<6p = 0.17). Times t are given in units of 
T= 10 4 r . 

to cylindrical and spherical structures. Consequently, the 
copolymers aggregate to rods and/or spheres in most of 
the parameter region. Close to the CMC, however, more 
complex structures are formed: Ring micelles, toroidal 
micelles, and even cage- like micelles. 

To understand why these complex micelles appear, one 
must inspect the pathways of structure formation in more 
detail. Both pathways to vesicle formation, PC and PG, 
are observed in our system (Fig.Q}. Likewise, rod forma- 
tion also proceeds via the two distinct pathways micelle 
coalescence (PC), or (anisotropic) micelle growth (PG), 
depending on the copolymer volume fraction (Fig. [2]): 
Coalescence takes place at $p > 0.2. Rings and toroidal 
micelles emerge at much lower copolymer volume frac- 
tion, and their formation is clearly driven by a growth 
mechanism. Fig. [SK shows the pathway to ring forma- 
tion, where spherical micelles first grow into small disks, 
a hole then nucleates at the center of the disks, and fi- 
nally, the perforated disks evolve into rings. Even closer 
to the CMC, the same mechanism leads to the formation 
of toroidal micelles: The micelles first grow into semivesi- 
cles or small vesicles, then several holes appear in the 
vesicle shells, until finally, the perforated vesicles grow 
into toroidal micelles (Fig [3b, c). The number of holes 
depends on the size of the embryo vesicle at the time of 
breakup. Close to the CMC, the dynamical stability of 
large vesicles increases: The initial number of micellar 
nuclei is small and they are far apart, hence the vesicles 
are not perturbed by the environment and break up late. 
As a result, large cage micelles can be obtained in the 
vicinity of the CMC (Fig. [3b). 

Next we examine the early stages of micelle aggrega- 
tion. To this end, we define the segregation parameter 
(7 = / dr\<frA(r)+$B(r)-<S> P \/V, where $aW and $p(r) 
are local monomer densities, and V is the volume. Look- 
ing at a as a function of time, we find that the copoly- 
mer segregation proceeds in an almost step-like fashion 
after a well-defined "incubation time" r* (Fig. [H inset). 
Fig.[4]shows the incubation times as a function of xbs for 
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FIG. 4: Inverse incubation time 1/r* as a function of xbs 
for different copolymer volume fractions x* as indicated. The 
values xbs have been shifted by x* = 1-3 — ll$p. The solid 
line corresponds to the function y = 5.3 • W~ e 'x 125 , and the 
dashed line shows y <x x^ 1 for comparison. Inset: Evolution 
of the order parameter a with time at xbs = 0, "^p = 0.17. 
The incubation time is defined as the time where the tangent 
of a(t) at the inflection point and the base line a = intersect 
each other. 



different copolymer volume fractions <&p. Remarkably, 
most curves collapse onto a single power law function 
1/r* = A ■ (xbs — X*) a with the exponent a w 1.25 after 
a simple shift of xbs by X*(^p) — 1-3 — ll$p. Only 
for the lowest copolymer volume fractions, $p = 0.13 
and $p = 0.15, do the data fail to collapse; however, 
the slopes of these curves, shifted by x* an d plotted in a 
double logarithmic way, are still comparable to a. 

A similar power law behavior has been observed pre- 
viously in a simulation study of vesicle formation in two 
dimensions [26] . It was explained in terms of the Cahn- 
Hilliard theory for spinodal decomposition, and x* was 
identified with the spinodal for macrophase separation 
between polymer and solvent. In the present case, how- 
ever, X*{^p) is significantly lower than that spinodal [32T|. 
it seems rather related to the CMC (cf. Fig. [2]). It is 
worth noting that data collapse and power law behavior 
is observed in a range of $p regardless of the final struc- 
ture and the dynamical pathway of structure formation. 
This suggests that the characteristics of the initial stage 
of segregation are universal and related to a spinodal- 
type instability. 

Once created, the micelles grow by attracting copoly- 
mers from the solution. In the pathway PC, the growth is 
supplemented by micelle coalescence. The question arises 
under which conditions this happens. In fact, most fu- 
sion events take place at early times (see Fig. Q]al,a2). 
At later stages, they are impeded by two factors: The 
formation of a well-segregated, swollen B-corona at the 
surface of the micelles, and the emerging copolymer de- 
pletion zone around the micelles. Hence the density of 
micelle nuclei at the early stage is a candidate quan- 
tity that might select between pathways. The crossover 
between pathways is observed at the copolymer volume 
fraction <f>p ~ 0.2. Taking into account that a large 
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fraction of copolymers assembles into droplets almost si- 
multaneously (at the time r*), and that these droplet 
nuclei roughly have the diameter 2R g , we can estimate 
the average distance D between droplets at given copoly- 
mer content^ via ( 2R g/ D ) 3 ~ $p- At $ P = 0.2, D is 
of the order R g . Hence the copolymers in solution are in 
contact with several droplets, they are attracted by all of 
them, and can serve as bridges that mediate fusion. 

In the last stage, the structures ripen. As long as they 
are still small (semivesicle state), classical Oswald coars- 
ening is observed (Fig. [T^3-a4, b3-b4): Small structures 
dissolve, large structures grow, driven by the competition 
of bulk and surface free energy. Once the structures have 
locally assumed their favored toroidal or bilayer struc- 
ture, the Oswald process stops and the ripening is gov- 
erned by much weaker thermodynamic forces, such as, 
e.g., those associated with the bending energy. The time 
scales of these processes are very slow and out of reach for 
our simulation method. Therefore, we have carried out 
a set of simulations using external potential dynamics 
with locally non-conserved (but globally conserved) den- 
sities. This dynamical model is less realistic, but much 
faster, such that we could also assess later stages of the 
aging process. Specifically, we studied the evolution of a 
system containing two vesicles with different initial sizes 
in two and in three dimensions. In the case of ring mi- 
celles or two dimensional vesicles, the contribution of the 
bending energy favors a uniform size distribution: The 
energy of a single ring of radius R is proportional to k/R 
(k being the bending rigidity), and the total energy for 
a fixed number of rings is minimal if all rings have the 
same diameter. Indeed, our two dimensional simulations 
showed that the sizes of the two rings converged in the 
course of the simulation: The small ring grew at the ex- 
pense of the large ring. For three dimensional vesicles, 
the situation is different: The bending energy of a vesicle 
is 47r/c, independent of its size. The total bending energy 
only depends on the number of vesicles, not on their size 
distribution. Consequently, we did not observe any sign 
of size uniformization, nor size disproportionation, in the 
three dimensional simulations. The same behavior has 
been observed in experiments [il]. 

To summarize, we have investigated the formation of 
toroidal micelles in copolymer solutions, and shown that 
such micelles may form in the vicinity of the CMC by 
a pathway that proceeds via the nucleation, growth, 
and subsequent breakup of vesicles. The competition 
of copolymer aggregation and self-assembly generates 
a kinetic trap that opens a route to manufacturing 
highly complex metastable structures. Several meth- 
ods to control these kinetic traps are conceivable: Con- 
trolling the number of initial nuclei by planting seeds 
(prenucleation) [iH HfJ, quenching from different initial 
states(e.g., quenching from a vesicle state), doping with 
additional components to tune the bilayer properties, or 



working with amphiphilic mixtures. 
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